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I. INTRODUCTION 



In this paper we generalize our previous study Jl| of chromatic polynomials and their 
asymptotic limits on strip graphs. The motivation for this work is that there are substances 
that exhibit nonzero ground state entropy, such as ice 0. A simple model that exhibits 
ground state entropy is the g-state Potts antiferromagnet (PAF) H|J for sufficiently large 
values of q. The ground state entropy is related here to the ground state degeneracy per site 
W by 5*o = kslnW. There is a direct connection with mathematical graph theory via the 
elementary equality Z(G,q,T = 0)paf = P(G,q), where Z(G,q,T = 0)paf is the partition 
function of the zero-temperature g-state Potts AF and P(G, q) is the chromatic polynomial 
PJH], giving the number of ways of coloring the n- vertex graph G with q colors such that no 
adjacent vertices have the same color, and the consequent equality 

W{{G},q) = lim P{G,q) 1/n (1.1) 

n—*oo 

where {G} denotes the n — > oo limit of the family of graphs G (e.g., a regular lattice, 
{G} = A with some specified boundary conditions). Since P(G,q) is a polynomial, there 
is a natural generalization of q from q G Z + to q G R+, and indeed to q G C. In general, 
W({G},q) is an analytic function in the q plane except along a certain continuous locus of 
points, which we denote B (and at possible isolated singularities which will not be important 
here). In the limit as n — > oo, the locus B forms by means of a coalescence of a subset of 
the zeros of P(G,q) (called chromatic zeros of G) M- ||T0f . 

In a series of papers [O]- two of us have calculated and analyzed W({G},q), both 



for physical values of q (via rigorous upper and lower bounds, large-g series calculations, 
and Monte Carlo measurements) and for the generalization to complex values of q. In Ref. 
U we presented exact calculations of chromatic polynomials and the asymptotic limiting W 
functions for strip graphs of a variety of regular lattices. For this purpose we developed and 
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applied a generating function method. Specifically, we considered strip graphs of regular 
lattices of type s of the form 

m 

(G s ) m j = ([[H)I (1.2) 

£=l 

where s stands for square, triangular, honeycomb, etc. Thus, such a graph is composed 
of m repetitions of a subgraph H attached to an initial subgraph / on one end, which, by 
convention, we take to be the right end. In the case of homogeneous strips, I = H. As is 
implicit in the above notation, we picture such strip graphs as having a length of L x vertices 
in the horizontal (x) direction and a width of L y vertices in the transverse (vertical, y) 
direction. Since for the limit of infinite length strips, the physical ground state degeneracy 
per site, W, and ground state entropy, So, are independent of the boundary conditions used 
in the longitudinal direction, we used free boundary conditions for technical simplicity in 



Ref. Jy. (Homeomorphic expansions of these strip graphs were also studied in Ref. [17].) 

In the present work we generalize our previous study by considering strip graphs of regular 
lattices in which there may be special subgraphs on both the left and right ends of the strip: 

m 

{G 8 )m,j,i = J([[H)I (1.3) 

1=1 

The end subsgraphs I and J may be different from each other and from the vertical rungs 
of the strip (thought of as a horizontal ladder). As before, we use the simplest type of 
boundary conditions - free ones - in the longitudinal direction, since in the limit of interest, 
namely that of infinite length, the physical ground state entropy is independent of these 
boundary conditions. However, since the size of the limit is finite in the transverse direction, 
it is of interest to consider both free and periodic boundary conditions in this transverse 
direction, and we do that here. Again, we use the generating function method. For the 
infinite-length limit of dl&Li H)I strip graphs, we proved that the nonanalytic locus B is 
independent of /, which is plausible, since in this limit, the subgraph I occupies a vanishingly 



small part of the full strip graph [|TJ. However, our study of the infinite-length limit of 
«/(n^=i H)I st r ip graphs reveals that in some cases, B does depend on the end subgraphs / 
and J. Furthermore, whereas for (fl^Li H)I strip graphs we found that B consists of arcs 
and possible line segments that do not enclose any regions in the q plane, in contrast, for 
J(n^Li H)I strip graphs, we show here that in certain cases B does enclose regions in the q 
plane. See also Ref. |TBJ and a companion paper on strip graphs of a particular Archimedean 
lattice |H§ denoted as (3 3 • 4 2 ). 



It will be convenient to define the function 



k-2 



where P(Ck, q) is the chromatic polynomial for the circuit (cyclic) graph Ck with k vertices, 

P(C k ,q) = (q-l) k + (-l) k (q-l) (1.5) 

This paper is organized as follows. In Section 2 we describe our present implementation 
of our generating function method. In sections 3-5 we apply this method to strip graphs 
of the square, honeycomb and triangular lattices. In section 6 we study strips of a certain 
Archimedean lattice known as the (4 • 8 2 ) lattice. Sections 7 and 8 present results for the 
square and triangular lattices with periodic boundary conditions in the transverse direction 
(equivalently, these may be called cylindrical boundary conditions). Some further discussion 
and conclusions are given in section 9. 



II. GENERATING FUNCTION METHOD AND END SUBGRAPHS 

We denote the chromatic polynomial describing the coloring of the strip graph (G s ) my jj 
with q colors as P((G s ) m jj,q). In the generating function method, this chromatic polyno- 
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mial is given by means of an expansion, in terms of an auxiliary variable x, of a generating 
function T(G S jj, q, x): 



T(G sM , q,x) = J2 P((G s ) m ,j,i, q)x m (2.1) 



m=0 

As before for strip graphs of the form ( |1.2|) , we find that the generating function T(G s ,j,i, q, x) 
is a rational function of the form 

r(G a ,j tI , q, x) = — — 2.2 

V(G S}JJ ,q,x) 

with (suppressing J, / dependence in the notation) 

Jmax 

M(G s ,q,x) = J2 a G s ,j(q)^ (2-3) 

3=0 

and 

kmax 

V(G s ,q,x) = l+J2b Gs Aq)x k (2-4) 

k=l 

where the ac a ,i and bc a ,i are polynomials in q. The degrees j max and k max of Af(G s , q, x) and 
T>(G S , q, x) as polynomials in x depend on the type and width of the strip G s . 

The method that we shall use for calculating the generating function is based on the 
addition-contraction theorem from graph theory. This was briefly described in Ref. [0, 
together with the equivalent deletion-contraction theorem; here we shall give a more detailed 
discussion. We first recall the statement of the addition-contraction theorem: let G be a 
graph, and let v and v' be two non-adjacent vertices in G. Form (i) the graph G a dd. in which 
one adds a bond connecting v and v', and (ii) the graph G contr . in which one identifies v and 
v'. Then the chromatic polynomial for G is equal to the sum of the chromatic polynomial for 
the graphs G a dd. and G contr .- For our discussion, we may begin by taking J and / to be the 
same as the other vertical rungs of the strip. By applying the addition-contraction theorem 
to the right-hand side of the initial H subgraph in (G s ) m , we obtain a set of strip graphs 



{G' B J) m with complete [pOj subgraphs labeled by j on the right-hand end. For example, let us 
consider a strip (G s ) m = YlT=i H wnere the width is such that to go from one transverse end 
to the other one traverses at least four vertices (including the two vertices on the transverse 
boundaries), and label the four vertices on the right-hand longitudinal end, in sequence, 
as V\, V2, t>3, U4. (For a strip of the square lattice with free transverse boundary conditions, 
this example corresponds to a strip of width L y = 4. For other types of strip this may 
correspond to different values of L y , as illustrated by explicit calculations below.) Now 
apply the addition-contraction theorem to the pair Vi,v^] this yields the equation 

P((G' s ) m -i, ?) = P((G' sM1A) ) m , ?) + P((G' s , vl= J m , q) (2.5) 

where {G' sb , x ^) m denotes the strip graph with the right-hand end modified by the addition 
of the bond connecting vertices v% and V4, and (G' sv v ) m denotes the strip graph with the 
right-end modified by identifying vertices V\ and v±. Applying the theorem again to each of 
these two strip graphs, one obtains the equation 

q)=P((G q)+P((G 
+P((G 

S,K2 (v 1 =V 3 ,V2 =Vi) ) m ' 

q)+P((G q) (2.6) 

Let us label the five resultant complete graphs on the right-hand end via the label i, with 
chromatic polynomials denoted by P((G' si )o, q). Thus, 

P((G' 8tl ) , q) = P(K,, q) = q(q - l)(q - 2)(q - 3) (2.7) 
P((G'J ,q)=P(K 3 ,q)=q(q-l)(q-2) for z = 2,3,5 (2.8) 

P({G' sA ) , q) = P{K 2 , q) = q(q - 1) (2.9) 

Next, working one's way leftward from the right-hand end, apply the addition-contraction 
theorem in a similar manner to the next set of transverse vertices, and label the five resul- 
tant complete graphs with the label j. Define a (square) matrix M (in this case 5x5) with 
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elements consisting of the chromatic polynomials of this set of graphs with complete sub- 
graphs i on the right and j on the left. This procedure transforms the initial strip into a sum 
of factorized strips with parts that overlap in complete graphs. Therefore the intersection 
theorem yields, for m > 1, 

P((G' s>j Uq) = E % S m "|' 9) = J2( MD hP((G' s ,) m -uq), (2.10) 
where D is a (5 x 5) diagonal matrix with elements 

D «=mk^) 1S!S5 (2 ' n) 

with the ordering as given above. (Do not confuse this matrix D with the denominator T> 
of the generating function T defined in eq. Q2.4j) . Because of (|2.6|) , the generating function 



of G' s can be written as 

T(G' s ,q,x) = J2 T (G' SJ ,q,x) (2.12) 

j 

Writing the generating function of G' s j in the form 

oo 

r(G' sJ ,q,x) = Y,P((G' s , J )m,q)x m - 1 (2.13) 

m=l 

and using ( |2.10|) with ( |2.13| ), we obtain 



r(G' s>j , g , x) = p((g' sj )i, g) + J2 x ( MZ} )^ E p (( G U)m-i, q)^ 2 

m=2 



= P((G , aJ ) 1 ,q) (MD) Jt T(G'^,q,x) (2.14) 

i 

The special case m = 1 in eq. (|2.10D yields 

P((G' S)i )i,g)=J] M ^ ( 2 - 15 ) 

i 

Let us define a vector V with elements T(G' S q, x), where the index j also labels the elements 
of the vector. Thus, eq. (|2.12j ) can be written as 



T(G' 8 ,q,x) = v b T -r (2.16) 

where v b is a vector with all elements equal to unity. In the example described above, 
Vb = (1, 1, 1, 1, 1) T . (Do not confuse the vectors Vb with the vertices Uj.) Similarly, let us 
define the vector P x with elements P((G' a j)x,q). Then eq. ( 2.15|) yields 



Pi = Mv a (2.17) 
where in this case v a = Vb- In vector form, eq. ( |2.14|) becomes 



T = V 1 +xMDT (2.18) 

Solving for T, we get 

r = (l-i¥D)" 1 Pi (2.19) 

Henceforth, for notational simplicity, we drop the prime on G and indicate the end subgraphs 
/ and J explicitly. Using ( [2.16| ) and ( |2.17| ) in eq. (|2.19| ), we obtain 



T(G sJ j, q, x) = v b T (l - xMD)- l Mw a = v b T x m (MD) m Mv a , (2.20) 



m=0 



where in this case, for the strip G s , v a = (1, 1, 1, 1, 1) T and v b = v a . Note that in general 
v a and Vb can be chosen independently, and different choices of these vectors correspond to 
different end graphs J and /, respectively. Thus, we can equivalently write (|2.1|) as 



P((G s ,j,i) m , q) = Vh T (M D) m Mv a (2.21) 
Factorizing the denominator of the generating function (|2.20|) , we have 

V(G a , JtI , q, x) = J](l - \ r (q)x), (2.22) 

r 

where the X r (q)'s are eigenvalues of the product of matrices MD. The subgraphs on the two 
ends of the strip (defined by v a and Vb) determine which eigenvalues enter in the product 



in equation (|2.22|) . Using this method, we have calculated chromatic polynomials and their 
asymptotic limits for various strip graphs of the form (G S) jj) m = J([\T=i H)I. We define 
(G s ,j,i) m =o as being equal to JI, where the end graphs / and J are, in general, different. 
For technical convenience, we also use a slightly different labelling convention than in Ref. 
|Ij], which can be illustrated as follows: for strips of the square lattice, for example, we take 
the end graphs I and J to refer to the respective sets of vertical bonds on the right and 
left ends of the strip, whereas in Ref. |TJ, / referred to the first vertical layer of squares. 
We now present some of these calculations and the analytic structure of the corresponding 
W({G sJJ },q). 

We start by defining the vectors which will be used to describe the subgraphs on the 
ends of the strips. Let us denote the four vertices on one end, in sequence, as V{, with 
i = 1, ...,4 and consider end subgraphs to be described by vectors v a and v b chosen from a 
set of possibilities that we list below: 

(i) 

Vl = (1,0,0,0,0) T (2.23) 

which corresponds to a strip with a complete graph K 4 on the end, as shown in Fig. |T|(a) 
for a strip of the square lattice of width L y = 4. This K& is formed by adding edges between 
vertices v% and V4, v\ and t>3, and v-i and v 3 . 

(ii) 

v 2 = (0,1,0,0,0) T (2.24) 

corresponding to a strip with vertices V\ and V3 contracted and an extra edge connecting 
vertices v 2 and w 4 , yielding a K 3 (vi = v 3 ) on the end. Fig. |l](b) illustrates this for a strip of 
the square lattice of width L y = 4. 
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Ill 



v 3 = (0,0, 1,0, Of (2.25) 

which corresponds to a strip with vertices v 2 and contracted and an extra edge connecting 
vertices v% and v^, yielding a K 3 (v 2 = 1)4) on the end. Fig. |l](c) illustrates this for a strip of 
the square lattice of width L y = 4. 



,iv 



v 4 = (0,0, 0,1, Of (2.26) 

which corresponds to a strip with vertices V\ contracted with w 3 and v 2 contracted with V4, 
forming a K 2 {y\ = v 3, v 2 = V4) on the end. Fig. |l](d) illustrates this for a strip of the square 
lattice of width L y = 4. 

(v) 

v 5 = (0,0,0,0,lf (2.27) 

which corresponds to a strip with vertices v\ and V4 contracted, yielding a K${v\ = V4) on 
the end. Fig. |I|(e) illustrates this for a strip of the square lattice of width L y = 4. 



vi 



v 6 = (l,l,0,0,0f (2.28) 

which corresponds to a strip with one extra edge connecting vertices V\ and V4 and another 
connecting vertices v 2 and t> 4 . Fig. |](f) illustrates this for a strip of the square lattice of 
width L y = 4. 

(vii) 

v 7 = (l,0,l,0,0f (2.29) 



which corresponds to a strip with one extra edge connecting vertices v\ and v 4 and another 
connecting vertices V\ and v^. Fig. |l](g) illustrates this for a strip of the square lattice of 
width L y = 4. Note that for a strip of the square lattice this end subgraph is equivalent 
to that described by vector V6 = (1, 1,0, 0,0) T . However, for other types of strips V7 may 
represent a different end subgraph. 

(viii) 

v = (1,1, 1,1, If (2.30) 

where the end subgraph is the same as the vertical rungs of the strip. Fig. |](h) illustrates 
this for a strip of the square lattice of width L y = 4. Fig. [I] illustrates these end subgraphs 
for a strip of the square lattice of width L y = 4. The different subgraphs are shown on the 
right end of the strip. Similar subgraphs can be considered on the left end of the strip. In 
the remainder of this section we consider the end subgraphs of the strips described by the 
set of vectors above, eqs. fl2.23| )-( fOU| ). 
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FIG. 1. Illustration of end subgraphs on a strip of the square lattice of width L y — 4 and free transverse 
boundary conditions. Cases (a), (b), (c) and (d) contain non-planar end subgraphs, while cases (e), (f), (g) 
and (h) contain planar end subgraphs. 

We label free and periodic boundary conditions in the transverse (y) direction as 

FBC y , PBC y (2.31) 

respectively, and similarly for the longitudinal (x) direction. The strips have FBC^ and, 
unless otherwise stated, also FBC y . For strip graphs with FBC X and FBCy, we shall restrict 
ourselves to cases where the product of matrices MD is a 5 x 5 matrix. For strip graphs with 
FBC X and PBC y we shall restrict ourselves to cases where the product of matrices MD is a 
4x4 matrix. These are strip graphs with cylindrical boundary conditions, with transverse 
cross sections comprised by four vertices and four edges, forming squares. 

Let us denote the four vertices belonging to the transverse cross section of the strip on 
one of the longitudinal ends of the strip by v iy with i — 1, ...,4, in sequence. The product of 
matrices MD is a 4 x 4 matrix in this case and vectors describing end subgraphs can be 

(i) 

Vl = (1,0,0, Of (2.32) 

which corresponds to a strip with a K4 on the end. This K4 is formed by adding edges 
connecting vertices v± and v% and vertices t>2 and V4. 

(ii) 

v 2 = (0,l,0,0) T (2.33) 

corresponding to a strip with vertices v 1 and t> 3 contracted and an extra edge connecting 
vertices vi and V4, yielding a K 3 (vi = 1*3) on the end. 

(iii) 
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v 3 = (0,0,1, Of (2.34) 

which corresponds to a strip with vertices v 2 and v& contracted and an extra edge connecting 
vertices V\ and V3, yielding a K 3 (v 2 = V4) on the end. 

(iv) 

v 4 = (0,0,0, 1) T (2.35) 

which corresponds to a strip with vertices v 2 and V4 contracted and vertices v\ and V3 
contracted forming a K 2 {v\ = v 3 ,v 2 = V4) on the end. 

(v) 

v 5 = (1,1,0, Of (2.36) 

where the ending subgraph is equal to the repeating unit with an extra edge connecting the 
non- adjacent pair of vertices v 2 and 1*4. 

(vi) 

v 6 = (1,0,1, Of (2.37) 

similar to (v), with the extra edge connecting the other pair of non-adjacent vertices, namely, 
V\ and i> 3 . 

(vii) 

v= (1,1,1, If (2.38) 

where the ending subgraph is equal to the vertical rungs along the strip. 

In our studies of strip graphs with PBC^ we only consider subgraphs on the two free ends 
of the strips described by the vectors listed in eqs. ( p.32|) -( ^73^ ). We next present results for 
a number of strip graphs of type G' s . 
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III. STRIP OF THE SQUARE LATTICE OF WIDTH L y = 4 

For a strip of the square (sq) lattice of width L y = 4 vertices and FBC y , we have (with 
f S q,i = fi to save space) 



( (g-3)/ /e (g-3)/ /i (g-3)/ /i (g-3)/ / 2 (g-3)/ /i 

(9 - 3)/ /i /o/s (?-2)/ / 2 / / 5 



M 



(?-3)/ / 2 
V(?-3)/o/i 



where 



(?-3)/„/i (g-2)/o/ 2 



/0/5 
/0/4 



/0/4 

/0/3 /0/5 /0/4 

/o/ B g(g-l)£>4 (g 2 -5g + 7)/ 
/0/4 (g 2 -5g + 7)/ /0/3 y 



(3.1) 



/s 9 ,o = 9(9- !)(<?- 2 ) 



(3.2) 



= <? 3 - 7g 2 + 19g - 20 



(3.3) 



/«g,2 = g 2 - 5g + 



/ a ,,3 = g 3 - 6g 2 + 14g - 13 



(3.4) 
(3.5) 



/ ag ,4 = g 3 - 7g 2 + 18g - 17 



(3.6) 



/ sg ,5 = g 2 - 4g + 5 



(3.7) 



/ s ,,6 = g 4 - 10g 3 + 41g 2 - 84g + 73 



(3. 



Hence, using eq. (|2.11|) with eqs. ( [2.7|) - (|2.9|) , we obtain the matrix product MD: 
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MD 



/ r 

h 


(q 


-3)/i 


(9 


-3)/i 


(g-2)(g-3)/ 2 


(?-3)/i ^ 


h 




fs 


(5 


-2)/ 2 


(g-2)/ 5 


/4 


h 


(9 


-2)/ 2 




/a 


(g-2)/ 5 


/4 


h 




/b 




h 




g 2 - 5g + 7 


U 




/4 




h 


(g-2)(g 2 -5g + 7) 


/a / 



(3.9) 



Three of the eigenvalues of this MD matrix, which we label as A,, % = 1,2,3 are the 
inverses of the three roots of the polynomial 



where 



6«,(4),i = ~q 4 + 7q 3 - 23g 2 + 41g - 33 
b sq (4),2 = 2g 6 - 23g 5 + 116g 4 - 329g 3 + 553g 2 - 517g + 207 



(3.10) 



(3.11) 
(3.12) 



?sg(4),3 



-q 8 + 16q 7 - 112g 6 + 449g 5 - 1130g 4 + 1829g 3 - 1858g 2 + 1084g - 279 (3.13) 



For all end subgraphs studied here, namely, considering any two pairs among fl2.23|) -( f2.30|) 
as the subgraphs on the two ends of the strip, these three eigenvalues contribute to the 
generating functions of the respective strips. The other two eigenvalues are A4 = 1 and 
A5 = (q 2 — 4g + 3). Depending on the end subgraphs chosen, either A4 or A5 or both 
contribute to the denominator of the generating function for the chromatic polynomial. 
However, we find that A4 and A5 are never leading eigenvalues. Hence, the nonanalytic locus 
B of W({G S q(4) t j t i},q) does not depend on the subgraphs J and I on the two ends of the 
strip and is given by Fig. 3(b) of Ref. [0. Of course, the actual chromatic polynomials and 
their zeros for specific finite strips do depend on J and /. 
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IV. STRIP OF THE HONEYCOMB LATTICE OF WIDTH L« = 3 



For a strip of the honeycomb (he) lattice of width L y = 3 and FBC y , we compute the 
MD matrix in the same manner as was discussed in detail above for the square lattice. We 
find that three of the eigenvalues of the MD matrix are the inverses of the roots of 

1 + &hc(3),is + hc(3),2X 2 + b hc{3)i3 x 3 (4.1) 

where the coefficients &^ c (3),i> &/ic(3),2 and bh c (3),3 are given by eqs. (3.14)-(3.16) in Ref. ([J; to 
save space, we do not list them here. For all end subgraphs studied here, namely, considering 
any two pairs among ( |2.23| )-( pllOD as the subgraphs on the two ends of the strip, these three 



eigenvalues contribute to the generating functions of the respective strips. The other two 
eigenvalues are identically zero. Therefore, just as was the case with the square lattice of 
width L y = 4, the denominator of the generating function, and thus £>, are not modified by 
the subgraphs on the two ends of the strip. 

V. STRIP OF THE TRIANGULAR LATTICE OF WIDTH L y = 4 

For a strip of the triangular (t) lattice with width L y = 4, we find that four of these five 
eigenvalues of the MD matrix are the inverses of the roots of the polynomial 

1 + b t (4),ix + b t{i):2 x 2 + b t{i):3 x 3 + & t ( 4 ),4^ 4 (5.1) 

where the coefficients &t(4),i, &t(4),2, &t(4),3 and &t(4),4 were given respectively by eqs. (B.15)- 
(B.18) in Ref. Since their expressions are somewhat lengthy, we do not list them here. 
Let us refer to these four eigenvalues as Aj, with i — 1, 2, 3, 4. For all end subgraphs studied 
here, namely, considering any two pairs among (|2.23|) -( p.30|) as the subgraphs on the two ends 



of the strip, these four eigenvalues contribute to the generating functions of the respective 
strips. The fifth eigenvalue is A5 = 1, which contributes to the generating function for certain 
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end subgraphs. A5 is leading in a region of the complex q plane and thus modify the curves 
B of non-analyticities of W({GW) j j}, g). In Table p] we show for various planar and non- 
planar end subgraphs J and /, whether or not A5 contributes to the generating function, 
and some features of the boundary B for each case. In turn, the presence or absence of A5 
in the generating function is determined by the end subgraphs I and J, or equivalently, the 
vectors v a and v& via the basic equation ( |2.21| ). Indeed, one sees from Table | that B encloses 
regions if and only if A5 appears in the generating function. This phenomenon, in which B 
can depends on the end subgraphs of the strip, is reminiscent of the dependence of certain 
aspects of complex-temperature phase boundaries and partition function zeros of the Potts 



model on boundary conditions found in Refs. p4|,[25 
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TABLE I. Illustrative end subgraphs and properties of generating function for a strip of the triangular 
lattice of width L y = 4. The notation J, I means that the vectors vj and Vi describe the end subgraphs on 
the two ends of the strip. P stands for planar and NP for non-planar. The third column shows whether or 
not A5 enters in the generating function. The fourth column lists some features of the boundary B. 



J, I 


1 *J- 

planarity 


does A5 enter: 


features 01 


V,V 


p,p 


N 


arcs 




P MP 


1\T 


arcs 


v, Vi, i = 5,6,7 


p,p 


N 


arcs 


Vi,v;, i = 1,2,3,4 


NP,NP 


Y 


arcs and one enclosed region 


Vi,v i; i = 5,6,7 


P,P 


N 


arcs 


vi,Vi, i = 2,3,4 


NP,NP 


Y 


arcs and one enclosed region 


vi, Vi, i = 5,6,7 


NP,P 


N 


arcs 


v 6 ,Vi, i = 1,2,3,4 


P,NP 


N 


arcs 


v 6 , v ; , i = 5,7 


P,P 


N 


arcs 
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In Fig. §(a) we show the analytic structure of the W function for an infinitely long strip 
of the triangular lattice with width L y = 4 and end subgraphs J, I that do not have A5 in 
the generating function. The resultant nonanalytic locus B is the same as that which we 
found in Ref. H (see Fig. 5(b) of that paper) for the infinitely long strip of the triangular 
lattice of the form (IX^Lj H)I with width L y = 4, for arbitrary /. The locus B in Fig. |2|(a) 
consists of arcs that do not separate the complex q plane into separate regions. This locus 
includes complex-conjugate multiple points of valence three; these are singular points on the 
algebraic curve constituted by B in the terminology of algebraic geometry [2~I|. Similarly, 



multiple points are observed in many of the loci B for strip graphs to be described below. 
Since the locus B forms by coalescence of zeros of the chromatic polynomial (called chromatic 
zeros of the graph) as the length of the strip graph goes to infinity, it is of interest to compare 
the locations of these zeros for a reasonably long finite strip with the asymptotic locus B. 
Accordingly, we show in this figure the chromatic zeros for a strip with I = H and length 
m = 12. 

In Fig. ^](b) we display the analytic structure of the W function for a strip of the 
triangular lattice with end subgraphs that yield A5 in the generating function. In this case the 
locus B is comprised by arcs which do enclose a self-conjugate region where A5 is the leading 
eigenvalue. The boundaries of this region cross the real q axis at q — qi — 3+ 2 v ^ = 2.618... 



and at q = 3. Following our earlier work |n,l2"|, we define q c as the maximal finite real value 
of q where W is nonanalytic for a given lattice, i.e., where B crosses the real axis. Hence, 
for the present strip, q c = 3. We note that for the infinite triangular lattice, q c = 4 [|i~2| . 
For comparison, in ||](b) we exhibit the zeros of the chromatic polynomial P((GfU\jj) m , q) 
with subgraphs on the right and left ends of the strip described by J = I = V2 and m — 10 
repeating units in the strip. 

We can understand this difference in the topology of the loci B as follows. One complex- 
conjugate pair of arc endpoints in the locus B in Fig. 0(a) occurs at q ~ 2.759 ± 0.154z. 



However, for the triangular L y = 4 strips with endgraphs that yield A5 in the generating 
function, this pair of would-be endpoints, and the contiguous portion of the associated arcs, 
do not actually occur on B because they lie inside the enclosed region in Fig. ||(b) where A 5 
is leading. Thus, as was anticipated in our discussion in Ref. |]], the difference in topology 
of B for the (II^i H)I strips and the present ^(11^=1 H)I strip of the triangular lattice 
is associated with the fact that for the former types of strips, the leading A's were always 
algebraic, whereas here a leading A, viz., A5 = 1 is a (constant) polynomial, rather than 
an algebraic function of q. Here we see how these features that B can (i) depend on J, / 
and (ii) enclose regions are both connected with each other via the basic eq. (12.21 ), which 
determines which A's appear in the generating function and, in particular, whether the 
leading A's are always algebraic, as for (fX^Lj H)I strips with nontrivial B loci, or can also 
include polynomial A's, as for ./(n^i H)I strips. 



19 




20 



FIG. 2. Analytic structure of the function W({G t (4),jj}, q) for a (a) (L x = oo) x (L y = 4) strip of the 
triangular lattice given by linim—,,^ J(Ilfci ^)I with end graphs /, J that do not yield A5 in the generating 
function (this includes the case where J is absent); (b) lim m _ +00 J(Il"li H)I with end subgraphs /, J that 
yield A5 in the generating function. For comparison, the zeros of the chromatic polynomials are shown for 
the strips with (a) I — H and m — 12 (whence n — 56); (b) subgraphs I — J described by v 2 and m = 10 
(whence n = 46). 

Thus, with changes in the subgraphs on the two ends / and J of a strip of the triangular 
lattice of width L y = 4, the analytic structure of W({Gt(i) t j,i}, q) changes from a set of arcs 
to structures which enclose regions. For this type of strip we only obtained enclosed regions 
when the subgraphs on both ends of the strip were non-planar. However, in Ref. JL9J an 
example is given where this happens with planar end subgraphs, albeit with a heteropolyg- 
onal Archimedean lattice, viz., the (3 3 • 4 2 ) lattice. Our present calculation shows that for 
the strip graphs of the form Jdl^i H)I, the nonanalytic locus B may depend on the end 
subgraphs / and J, in contrast to the property proved in Ref. Jl| that for strip graphs of the 
form dl^i the resultant B is independent of the end subgraph /. 



VI. STRIP OF THE (4 8 2 ) LATTICE OF WIDTH L y = 3 

We have also studied strips of lattices that are tilings involving more than one type of 
regular polygon. We recall p2]j23|JT6| JT[1 that an Archimedean lattice is a uniform tiling of the 
plane by one or more regular polygons, such that each vertex is equivalent to every other 
vertex. An Archimedean lattice is identified by the symbol 

i 

where the product is the ordered sequence of polygons that one traverses in making a com- 
plete circuit around a vertex in a given (say counterclockwise) direction, and the notation 
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indicates that the regular p-sided polygon pi occurs contiguously times in this circuit. 
In this section we consider a strip of the (4 • 8 2 ) lattice of width L y = 3 with free boundary 
conditions in the y direction. For this strip, the MD matrix is 5 x 5, and three of the 
eigenvalues of this matrix are the inverses of the roots of l + &488(3),i x + ^488(3), 2^ 2 + b^8(3),3X 3 , 
where the coefficients 6488(3) j, 3 — 1, 2, 3 were given in Ref. For all of the end graphs / 
and J studied here, namely considering any two pairs among (i)-(viii) in eqs. ( |2.23| )-( p.30| ), 
these three eigenvalues contribute to the generating functions of the respective strips. The 
other two eigenvalues are identically zero. Therefore, as was the case with the strip of the 
honeycomb lattice with L y = 3, the denominator of the generating function and thus B, are 
not modified by the end graphs / and J on the strip. 



VII. STRIP OF THE SQUARE LATTICE WITH PBC y 

In this section we consider a strip of square lattice with PBC^, i.e., periodic boundary 
conditions in the transverse direction (= cylindrical boundary conditions). We begin by 
observing that if one takes a strip of a regular lattice with FBC^ and identifies the vertices 
with lowest and highest values of y for each x to produce PBC y , the changes in B depend 
on the specific type of lattice and value of L y . For example, let us consider the square 
lattice with FBC y and L y = 4 and carry out the above-mentioned operation, equivalent to 
gluing together the upper and lower edges of the strip, so that the transverse cross sections 
are triangles. Although the infinite-length limit of the original strip with FBC^ yielded 
a nontrivial B, shown in Fig. 3(b) of Ref. Q, after the gluing operation, the chromatic 
polynomial factorizes trivially for the square strip with PBC^ and triangular cross section, 
and B — 0, i.e. there is no locus of points where W is nonanalytic in the q plane. In 
contrast, for a number of cases to be discussed below, the gluing operation takes an infinite 
strip with nontrivial B to another infinite strip with nontrivial and different B. As we shall 
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show below, for strips of the square (triangular) lattice with transverse cross sections forming 
squares, and end graphs / and J that are the same as the other vertical rungs of the strip, 
the resultant B does not (does) enclose regions, respectively. 

We begin by studying the strip of the square lattice with PBC y and transverse cross 
sections forming squares. Depending on one's labelling conventions, this corresponds to 
L y = 4 or L y = 5, where in the latter case, one interprets the periodic boundary conditions 
as identifying the top and bottom vertices for each value of x. It is interesting to observe 
that this graph can be regarded as part of a three dimensional simple cubic lattice. For end 
subgraphs I and J described by the vector v we obtain a generating function with j max = 1, 
k m ax = 2, and coefficients 

a sq (4), c ,o = q(q - 1)(<? 6 - H<? 5 + 55g 4 - 159g 3 + 282g 2 - 290g + 133) (7.1) 
a«K4),c,i = ~q(q ~ 1)(<? 2 - 3g + 3)(g 6 - 12g 5 + 61g 4 - 169g 3 + 269g 2 - 231g + 85) (7.2) 
6. g (4),c,i = ~q A + 8g 3 - 29g 2 + 55g - 46 (7.3) 

& sg (4), c ,2 = q 6 - 12g 5 + 61g 4 - 169g 3 + 269g 2 - 231g + 85 (7.4) 

In this case the MD matrix is the upper left-hand 4x4 submatrix of MD given above in 
eq. (|3.9p for the square strip with FBC r Two of the four eigenvalues of the MD matrix are 
the inverses of the roots of 1 + 6 sg (4) jC) ix + b sq ^)^2X 2 . These two eigenvalues, denoted Ai and 
A2, contribute to the generating functions for all end graphs J and / that we consider. The 
other two eigenvalues are A3 = 1 and A4 = (q — l)(q — 3), which contribute to the generating 
functions for some end graphs. In Table [ll| we show for various end subgraphs J, I whether 
or not A 3 and A 4 contribute to the generating function and some features of the boundaries 
B for each case. 
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TABLE II. Illustrative end subgraphs and properties of generating function for a strip of the square 
lattice with periodic boundary conditions on the transverse direction and transverse cross sections forming 
squares. Notation is as in Table Q, The second and third columns indicate whether or not A3 and A4, 
respectively, contribute to the generating function. 



J, I 




does A3 enter: 


does A4 enter; 
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V, V 




N 


N 






arcs 






Vi, Vi, i = 1, 


4 


Y 


N 
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Vi, Vi, i = 2, 


3 


Y 


Y 


arcs 


and 


one pair of enclosed 


rej 


dons 


Vi/vi, % = 5, 


6 


N 


Y 


arcs 


and 


one pair of enclosed 


rej 


dons 


v, Vi, i = 1, . 


,6 


N 


N 






arcs 






vi, v ; , % = 2,3,4 


v 
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6 


N 
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v 2 , v 3 




Y 


Y 
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rej 
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N 
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3 


Y 


N 
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6 


N 
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arcs 
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Y 


arcs 


and 


one pair of enclosed 




dons 
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FIG. 3. Analytic structure of the function W for an infinitely long strip of square lattice with periodic 
boundary condition on the transverse direction and transverse cross sections forming squares: (a) corresponds 
to end graphs I and J for which A4 does not enter in the generating function, (b) corresponds to cases where 
A4 enters in the generating function. For comparison, the zeros of the chromatic polynomial for a strip with 
m = 8 repeating units (n = 40 vertices in (a) and n = 36 vertices in (b)) and end graphs / and J described 
by (a) v [(b) v 2 ] on the right and left ends of the strip, are shown. 

We note that although A3 enters in the generating functions for some end graphs, it is 
never leading, so the analytic structure is not modified by its presence in the generating 
function. For end graphs where A4 is not present in the generating function the nonanalytic 
locus B of the W function, shown in Fig. 0(a), includes one complex-conjugate pair of arcs 
and one self-conjugate arc which crossing the real q axis at q ~ 2.3026. There is also a line 
segment on the real q axis from q = 2.2534 to q = 2.3517, which latter point is thus the 
value of q c for (the m —>■ 00 limit of) this strip. The locus B in Fig. 3(a) is somewhat similar 
to that which we presented in Fig. 3(a) of Ref. fl| for the infinitely long strip of the square 
lattice of width L y = 3 with FBC^ except that in that case, the self-conjugate arc crossed 
the real q axis at q — 2, which was thus the value of q c in that case. The value of q c that 
we find for PBC^, viz., q c ~ 2.35, is comparable to the value q c ~ 2.27 that we found for 
the infinitely long strip of the square lattice with FBC y (see Fig. 3(b) of Ref. I]) and width 
L y = 4. The endpoints of the arcs and the line segment correspond to the branch points of 
the eigenvalues Ai and A2. For end graphs I and J where A4 is present in the generating 
function the analytic structure of the W function, shown in Fig. 0(b), has an extra pair 
of complex-conjugate regions. In these enclosed regions the eigenvalue A4 is leading. The 
complex-conjugate pair of branch points of the eigenvalues Ai and A2 at q = 1.992 ± 1.594i 
lie within the complex-conjugate enclosed regions, where A4 is leading. Hence, these branch 
points do not correspond to arc endpoints. 
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VIII. STRIP OF THE TRIANGULAR LATTICE WITH PBC 



We next consider a strip of the triangular lattice with PBC^ and transverse cross sections 
forming squares. As above, depending on one's labelling conventions, this corresponds to 
L y = 4 or L y = 5, where in the latter case one identifies the top and bottom vertices for each 
value of x. The subgraphs on the horizontal ends of the strip are described by the vector v. 
We obtain a generating function with j max = 1, k max = 2, and coefficients 

Ot(4),c,o = q(q - 1)(? - 2)(g - 3)(g 4 - 10g 3 + Alq 2 - 84g + 71) (8.1) 
a t( 4) )C ,i = -2q(q - l)(q - 2){q - 3) 3 (g 2 - 5q + 5)(g 2 - 3q + 3) (8.2) 
6t(4) )C) i = ~(q - 3)(g 3 - 9q 2 + 33q - 48) (8.3) 

h m , c ,2 = 2(q -2){q- 3) 3 (g 2 - 5q + 5) (8.4) 

Two of the four eigenvalues of the MD matrix are the inverses of the roots of 1 + 6t(4), c ,ix + 
6 4 ( 4 ) iCj2 x 2 . These two eigenvalues, denoted Ai and A 2 for the + and — signs of the square 
root, respectively, contribute to the generating functions for all end graphs I and J that we 
consider. The other two eigenvalues are A 3 = and A 4 = 2, the latter contributes to the 
generating functions for some end graphs. In Table [TT1] we show for various subgraphs J, I 
on the two ends of the strips, whether or not A4 contributes to the generating function and 
some features of the boundary B for each case. 
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TABLE III. Illustrative end subgraphs and properties of generating function for a strip of the triangular 
lattice with periodic boundary conditions on the transverse direction and transverse cross sections forming 
squares. The notation is as in Table |[ The second column indicates whether or not A4 contributes to the 
generating function. 
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FIG. 4. Analytic structure of the function W for the strip of triangular lattice with periodic boundary 
condition on the transverse direction, transverse cross sections forming squares, and length L x — oo. (a) 
corresponds to end graphs / and J for which A4 does not enter in the generating function, (b) corresponds 
to cases where A4 enters in the generating function. For comparison, the zeros of the chromatic polynomial 
for a strip with m = 8 repeating units (n = 40 vertices) and end graphs / and J described by (a) v [(b) Vi] 
on the right and left ends of the strip, are shown. 

For end graphs where A4 is not present in the generating function the analytic structure 
of the W function, shown in Fig. f|(a), is formed by one pair of complex-conjugate arcs 
and one self-conjugate region. Outside (inside) the enclosed region X 1 (A 2 ) is leading, The 
endpoints of the arcs correspond to the branch points of the eigenvalues Ai and A 2 . The 
boundary B crosses the real q axis at q ~ 3.481 and at q — 4. For end graphs where A4 is 
present in the generating function the analytic structure of the W function, shown in Fig. 
f|(b), is qualitatively the same as the previous case, except that the enclosed region is larger 
and A4 is leading therein. The boundary B crosses the real q axis at q = 3+ 2 v ^ = 2.618... and 
at q = 4. It is interesting to note that for a strip of the triangular lattice with PBC y and 
transverse cross sections forming squares and different end graphs / and J, our results yield 
a value of q c which coincides with the value for an infinite triangular lattice, namely, q c = 4 



T4 . However, this is not a general result; in our analysis above of the strip of the square 
lattice with PBCy and transverse cross sections forming squares, we obtained q c ~ 2.35, 
which is somewhat less than the value q c = 3 for the infinite square lattice. 



IX. DISCUSSION AND CONCLUSIONS 

In this paper we have presented exact calculations of chromatic polynomials P and asymp- 
totic limiting W functions for strip graphs of the form J(n^i where J and / are various 
subgraphs on the left and right ends of the strip, and (n^Li ^0 are strips of various regular 
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lattices consisting of m-fold repetitions of subgraph units H. We have also studied the ef- 
fects of using periodic as well as free boundary conditions in the transverse direction. For 
the respective strip graphs we have determined the loci B where W is nonanalytic in the 
complex q plane. Our present analysis generalizes our earlier calculations for strip graphs of 
the form ([\f =1 H)I. 

In contrast to the dl^i H)I strip graphs, where B (i) is independent of /, and (ii) 
consists of arcs and possible line segments that do not enclose any regions in the q plane, 
we find that for some J(f7^ 1 H)I strip graphs, B (i) does depend on / and J, and (ii) can 
enclose regions in the q plane. We have explained these differences in the present work and 
have related them, via eq. ( p. 21 ) to whether the leading A's are algebraic, as in the case of 



(n^i H)I strips with nontrivial B, or can also include polynomial A's, as in certain of the 
J(n^i H)I strips. It should be noted that even if all of the leading A's are algebraic, the 
resultant locus B can still enclose regions. This can occur when this locus B is noncompact 



in the q plane, extending infinitely far from the origin [14]. Indeed, we showed in Ref. [14| 
that whenever B is noncompact in the q plane, passing through 1/q = 0, it encloses regions. 
An example of a family of graphs with only algebraic A's which has a locus B that encloses 
regions (and is noncompact) is the {Uk} family discussed in detail in (the third paper of) 
Ref. |T3| . Our present study thus elucidates the effects of different end subgraphs I and J 



and of free versus periodic transverse boundary conditions on the function W describing the 
ground state degeneracy, per site, of the Potts antiferromagnet on infinitely long 2D strips 
of various lattices. 
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